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Abstract 

Any account of the emergence of classicality from quantum theory must address the fact that the 
quantum operators representing positions and momenta do not commute, whereas their classical 
counterparts suffer no such restrictions. To address this, we revive an old idea of von Neumann, and 
seek a pair of commuting operators X, P which are, in a specific sense, "close" to the canonical non- 
commuting position and momentum operators, x,p. The construction of such operators is related 
to the problem of finding complete sets of orthonormal phase space localized states, a problem 
severely limited by the Balian-Low theorem. Here these limitations are avoided by restricting 
attention to situations in which the density matrix is reasonably decohered (i.e., spread out in 
phase space). 

PACS numbers: 
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I. INTRODUCTION 



A. Opening Remarks 



The canonical commutation relation between position and momentum operators, 



[x,p] = ih 



(1.1) 



lies at the heart of all quantum phenomena. Physically, this relation corresponds to the fact 
that measurements of position and momentum depend on the order of the measurements. 
It is this relation, in essence, that is responsible for the key differences between classical and 
quantum mechanics, since in classical mechanics, measurements of position and momentum 
can be made in such a way that the order makes no difference. It follows that any account 
of the emergence of classical behaviour from quantum theory must reconcile these two very 
different aspects of classical and quantum theory. A typical point of view is that classical 
mechanics emerges only at a very coarse-grained level and for sufficiently coarse-grained 
samplings of position and momentum, their non-commutativity makes little difference [l|. 

Still, one wonders whether there is a deeper or more precise way of reconciling the non- 
commuting quantum operators with their commuting classical counterparts. Indeed, this 
clearly troubled the founders of quantum theory, since von Neumann addressed the issue in 
his 1932 book. Mathematical Foundations of Quantum Mechanics {2!. He noted that when 
we make observations of a macroscopic system, we are in fact able to make observations of 
position and momentum simultaneously (although imprecisely, of course). This suggested 
to him that the measurements we make do not in fact correspond directly to the usual 
operators x and p, but to some other operators, X and P, say, which commute 



and which must in some sense be "close" to the original operators £, p. Such a pair of 
operators could be particularly useful in bringing a degree of precision to discussions of 
emergent classicality. The aim of this paper is to discuss the construction and utility of such 
operators. 



[X,P]=0 



(1.2) 
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B. The Imprecision of Classical Physics 

To what degree do macroscropic measurements fix the position and momentum opera- 
tors in quantum mechanics? When we make macroscopic measurements of, say, a particle, 
there will be imprecisions Ax, Ap in the specifications of position and momentum. These 
imprecisions may be "small" to classical eyes but they will typically be very large compared 
to the quantum scale. It will be useful for what follows to get a quantitative idea of this. 
Suppose that the position is measured to a precision of 10~^m and the velocity to within 
10~^m/s (both extremely precise specifications, by macroscopic standards). For a mass of, 
say, 10~^kg, we then have, 

^ ~ 10- (1.3) 
n 

This means that when we specify the phase space location of a classical system in a way that 
is very precise to classical eyes, at the quantum scale it could be in any one of about 10^^ 
phase space cells. There is therefore a considerable amount of freedom at the quantum level 
to redefine the position and momentum operators without making any noticeable difference 
to macroscopic observations. Perhaps within this freedom there is the possibility to find 
position and momentum operators which commute. 



C. Motivations 



What are the motivations for constructing commuting operators which are close to the 
position and momentum operators? An important indicator of emergent classicalityis ap- 
proximate diagonality of the density operator p (often referred to as decoherence) [3!]. For 
example, for a point particle interacting with a thermal environment, it may be shown that 
after a short time scale, the density operator approaches the form 



p = J dpdq f{p,q) \iJpg){^pg\ (1.4) 
where f{p,q) is a non-negative function and {ippq) are a set of phase space localized states 



(such as generalized coherent states) ^|. This indicates that it is approximately diagonal 
in both position or momentum. (It cannot be exactly diagonal in both, since they do 
not commute.) This means that there is negligible interference between different values 
of position and momenta and, loosely speaking, they may then be treated as if they are 



3 



classical. The statement is imprecise, however, since there is still some interference, so the 
variables are only imprecisely defined. 

One would like to be able to make a more exact statement about diagonality of p and this 
is possible in terms of the commuting operators X, P. For example, using the eigenstates 
\ipnm) of X, P one could construct a density operator of the form 

P = y^Jnmli^nm) {'ipnml (1-5) 
nm 

which is exactly diagonal in X and P, indicating there is exactly zero interference between 
different values of these quantities. One can then choose the coefficients fnm to make p' as 
close as possible to p. The question of the degree of approximate diagonality is therefore 
shifted to the question of the closeness of the density operator to a pseudo-classical density 
operator p'. 

A more comprehensive approach to emergent classicality is the decoherent histories ap- 



proach 



10| . There, the central object of interest is the decoherence functional. 



= Tr {PaM ■ --PaMpPa'Stl) " ' ' Pa'^iQ) (1-6) 



The histories are characterized by the initial state p and by the strings of projection operators 
Pa(t) (in the Heisenberg picture) at times ti to t„ (and a denotes the string of alternatives 
«! ■ ■ ■«„). Intuitively, the decoherence functional is a measure of the interference between 
pairs of histories a, a'. When 

ReD(a,a')=0 (1.7) 

for a ^ a', we say that the histories are consistent and probabilities p{a) = D(a,a) obeying 
the usual probability sum rules may be assigned to them. (Typically, the physical mecha- 
nisms producing consistency actually cause the stronger condition -D(a, a') = for a ^ a' 
to be satisfied, which is referred to as decoherence.) One can then ask whether these prob- 
abilities are strongly peaked about trajectories obeying classical equations of motion, and if 
they are, we can say that the system is emergently classical. 

For histories in which the projections Pa^itk) are onto positions at different times of 
a point particle interacting with a thermal environment, the decoherence functional, like 
the density operator, is approximately diagonal [g], [io|. The approximation is typically 
exceptionally good, but still only approximate. Again one wonders whether more exact 
statements can be made. Indeed, it has been conjectured that approximately consistent 



histories can be in some sense distorted into exactly consistent ones [ll|] . There are a number 
of ways in which a set of histories could be distorted: one can change the initial state, the 
times of the projections, the widths of the projections or the operators who spectrum is 
projected onto. The possible existence of the commuting variables P, X suggests a particular 
way of distorting the histories so as to make them exactly decoherent. 

Each position at time t is, in the Heisenberg picture, a function of the canonical pair, 
so Xt = ft{x,p). The projections onto Xt at different times do not commute, since 
7^ in general. The decoherence functional is therefore not diagonal in general, but 
can be approximately diagonal if the system is coupled to a thermal environment. 

Now suppose we consider projections onto the variables Xt = ft{^, P) at different times. 
Under reasonable dynamics, Xt will be close to Xt as long as x,p are close to X,P. The 
operators Xt do commute at different times so all the projections commute and, as is easy 
to see, the decoherence functional will be exactly diagonal. So exact decoherence is achieved 
by shifting the operators x,p to the commuting pair X,P. Furthermore, it is known that 
the probabilities for histories of positions Xt are typically peaked about classical evolution. 
The probabilities for the commuting variables Xt will therefore have the same property if 
the Xt are close to Xt- The question of approximate decoherence and emergent classicality 
is therefore shifted to the question of the closeness of the old and new operators. 

These then, at least in outline, are the reasons why a commuting pair of position and 
momentum-like operators may be useful for discussing emergent classicality. We turn now 
to the construction of such operators. 

D. Von Neumann's Construction 

Von Neumann outlined a prescription whereby the commuting operators X and P may 
be constructed [2^. This involved first taking a discrete subset \m,n) of the coherent states 
\p,q), with one state per cell of size 2'7ih (a von Neumann lattice). He alleged that these 
states are complete (this was later proved [l2, [l3, Q). He then stated that they may be 
orthogonalized using the Schmidt process to produce an orthonormal set \ipnm)- From these, 
he constructed position and momentum-like operators 

X = ^na\i)nm){'ipnm\, -P = ^ | ^nm) (^nm | (1-8) 
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(where a is a constant with the dimension of length). These operators clearly commute. He 
argued that these new operators are indeed "close" to the old ones, in the sense that 

(V'nml(x-X)Vnm) (^nm | (p - P) Vnm) < K^^^ (1-9) 

Von Neumann's calculations imply that the constant K is about 1, 800 (but he thought more 
detailed calculations could give a smaller value). 

However, von Neumann's prescription is at best at sketch of how this works and he 
certainly did not give full details (such as the explicit form of the lipnm))- Furthermore, we 
now know a lot more about phase space localized states than was known in 1932, and, as will 
be described below, there are obstructions to constructing such states. These obstructions 
do not necessarily apply to what von Neumann did, but nevertheless, it is still interesting 
to revisit his ideas from a more modern perspective. 



E. A General Approach and the BaUan-Low Theorem 

We start with a more general statement of the problem. We consider a set of states of 
the form 

l^nm) = Unmli^) (1.10) 

where {ip) is a fiducial state and Unm is the unitary shift operator, 

Unm = exp ( —nap — —mbx ) (l-H) 



h h 

A particularly interesting case is that of a von Neumann lattice, in which case h = 2TTh/a 
and the translations in the p and x directions then commute, and we have 

Unm = (-1)"" exp (^^nap^ exp (^-t'^x^ (1.12) 

There is then one state per cell of size 27rh, as in von Neumann's case. It is of interest to find 
a fiducial state in Eq. fll.lOp such that the states are complete and orthonormal, that is, 

\^nm){^nm\ = 1 (1-13) 

nm 

{■ipnml^n'm') = ^nn'^mm' (1-14) 

Fiducial states leading to states satisfying these properties are easily found and we will 
exhibit a set below. There is, however, a crucial difficulty. According to the theorem of 



Balian and Low 
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17|, if the three properties (11. 101) . (11.131) and (I1.14p are satisfied, 



then the fiducial state has either (Ax)^ or (Ap)^ infinite, so is not phase space locahzed. 
If we used such states to construct commuting operators X,P as von Neumann did, then 
at least one of the averages in Eq. (11.91) would diverge, and there would be no sense in which 
the new operators are close to the old. (Note that von Neumann claims to have used the 
Schmidt procedure to construct his orthonormal set, which one would not expect to produce 
states satisfying Eq. Ol.lOp . so his construction does not necessarily fall foul of the Balian-Low 
theorem). 

The problem of constructing orthonormal phase space localized states is one of great 
interest in a number of fields so some effort has been expended in finding ways around the 
Balian-Low theorem. Zak has proved some interesting results in this area. In Ref |18j], he 
showed that the coherent states restricted to a von Neumann lattice |m, n) obey a sort of 
orthogonality relation if the usual inner product (m, n|m', n') is averaged over a single phase 
space cell. It is not yet clear if this result can be used to produce commuting position and 
momentum operators. He has also considered complete orthonormal sets of states which 
are localized in position, but double-peaked in momentum (so localized in p^, but not in p. 
j'rom these one can construct commuting operators X and P^, which are "close" to x and 
lol . This is tantalizing close to the goal of this paper, but not quite there (and also suggests 
that the p —p transformation plays a crucial role in the Balian-Low theorem). Many of 
these and similar results are proved using the so-called kq representation, a technique which 
is particularly well-adapted to these problems [20|]. 

To avoid the Balian-Low theorem, one has to drop one of the three requirements (11.101) . 
(I1.13P and (11.140 in order to get phase space localization. For the purposes of this paper, 
which is to construct useful commuting position and momentum operators, the orthogonality 
Eq. (ll.l4p is essential. We will therefore explore the possibility of dropping the other two 
requirements. Eq. Ol.lOp . the requirement that the states be obtained by translation of a 
single fiducial state is mainly for practical convenience, so there is no harm in relaxing this 
as long as the resulting states are not unmanageable. 

More significantly, we will relax the requirement of completeness, Eq. (ll.l3l) . The motive 
behind this is as follows. We are interested in using commuting operators to discuss emergent 
classicality. In practice, this means that we are only concerned with the physical situation 
in which the density matrix of the system has undergone a degree of decoherence. This 
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means that it is approximately diagonal in both position and momentum, or equivalently, 
its Wigner function is reasonably spread out in phase space. The density matrix is therefore 
insensitive to the fine structure of phase space and it seems reasonable to suppose that the 
physics could be well-described by a less than complete set of states, if they are carefully 
chosen. Of course, this is a quantitative matter that needs to be checked in detail and we 
will do this. 

For the purposes of constructing commuting position and momentum operators, we will 
actually use a construction slightly more general than the one indicated by von Neumann 
in Eq. fll.Sp . In particular, we will look for commuting operators X, P of the form 

X = J^XnEnm (1.15) 

nm 

P = Y^PmEnm (1-16) 

nm 

Here, X„ and Pm are c-numbers and Enm are projection operators localized onto a region of 
phase space labelled by ra, m which will consist of more than one 27r^-sized cell. They are 
exclusive 

EfimEn'm' E^m^nn'^mm' (I'l"^) 

and exhaustive 

J2Enm = l (1.18) 

nm 

We will also insist that they are obtained from unitary shifts of a single projector E, 

Enm UfimE iUnni) (1.19) 

Here, Unm shifts from one cell to the next, so is of the form Eq. fll.llj) . with a, h chosen so that 
the translation in the position and momentum directions commute, but ah > 2TTh. These 
three conditions are the natural generalization for projection operators of the requirements 
( ll.lOp . (11.131) and (11.141) . and the original case is obtained with the choice 

Enm = |^nm)(^nm| (1-20) 

The quantity TtE is a measure of the number of 27r^-sized cells projected onto, so clearly 
TtE = 1 in the pure state case, but TiE ^ 1 in our case, as we will see. 

One might have thought that this more general construction could avoid the Balian-Low 
theorem, since the restrictions on Enm are in fact weaker than the restrictions (ll.lOp . (I1.13p 
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and fll.l4p for pure states. This is not in fact the case. We will prove a modest extension 
of the Balian-Low theorem which shows that there is no phase space localized projector 
E satisfying the three requirements Eqs. fll.17l) . fll.181) and (I1.19p . However, what we will 
do is find an "almost" complete set of phase space localized pure states from which we 
can construct a projector E that satisfies the exhaustivity condition Eq. fll.lSp to a good 
approximation when acting on sufficiently decohered density operators. From this we can 
construct commuting operators X and P may with useful properties. 



F. Earlier Work 



Finally, we briefly mention a related approach. An earlier 



attempt to construct commut- 



ing position and momentum operators was considered in Ref. |2l|. This construction involved 
doubling the orig inal Hilbert space and then using operators deflned on this enlarged space. 
(See also Ref. 22| for similar ideas). The resulting theory is essentially the same as 't Hooft's 
deterministic quantum theory 23|. However, this is no longer standard quantum theory. In 
the present work, by contrast, we stay within the framework of standard quantum theory. 



G. This Paper 

In Section 2, we briefly summarize some known properties of the Wigner function which 
help to make precise the idea that the state is sufficiently spread out in phase space. We also 
briefly note that the Wigner function naturally suggests an alternative method of deflning 
commuting position and momentum operators. In Section 3, we prove a modest general- 
ization of the Balian-Low theorem for projector operators. In Section 4 we introduce an 
orthornormal set of phase space localized states that are "almost" complete. We then use 
them to construct a set of phase space localized projection operators which are almost 
exhaustive. In Section 5 it is shown that the incompleteness does not matter if the density 
operator of the system is reasonably spread out in phase space. In Section 6, we use the 
projection operators Enm to construct commuting position and momentum operators and 
compute the "distance" between these operators and the usual canonical pair, as in Eq. (11.91) . 
In Section 7, we compare the probabilities for X and P with those for the usual position and 
momentum operators and flnd them to be close. We summarize and conclude in Section 8. 
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This contribution is based on the pubhshed work Ref. 



II. SOME PROPERTIES OF THE WIGNER FUNCTION 

It will be useful for the rest of the paper to briefly summarize here some aspects of 
quantum mechanics in phase space and the Wigner function. Most of this is standard 
material and may be skipped by the informed reader, except the brief comments at the end 
of this section. 

The Wigner representation for a density operator p (or indeed a wide class of operators) 
is defined by 

W{p, <l) = l^j di e-i^^ p{q + \i.q- \i) (2.1) 

with inverse 

p(x, y) = j dp ei^^-y^ Wip, ^) (2.2) 

Many calculations involving operators are usefully expressed in the Wigner representation. 
For example, 

Tr (AB) = 2'Kh j dpdq Wa{p, q)WB{p, q) (2.3) 

where Wa{p, q) and Wb{p, q) are the Wigner functions of A and B. 

We will be interested in later sections in the behaviour of the density operator or Wigner 
function in a simple model of the decoherence process. We take the simplest case of a single 
particle coupled to a thermal environment in the limit of high temperature and negligible 
dissipation, with no external potential. The master equation for the density matrix p(x, y) 
is, 

dp ih ( (9V D ^ 



dt 2m \dx'^ dy"^ J h? 
where D = 2m'jkT. In the Wigner representation, the corresponding Wigner function obeys 

the equation, 

dVV _ _p_dVV ^ 
dt m dx dp^ 
Evolution according to the master equation Eq. (12.41) tends to produce approximate diago- 

nality in position and momentum. In the Wigner representation, this appears as a spreading 

out phase space. Indeed, using Eq. (12.51) it may be shown that 

(Ap)2 = 2Dt+{/\p)l (2.6) 

2 Dt^ t^ 2 

iMt = + (^P)l—2 + + (^^)o (2-7) 
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h ^ \ h ) 



where 

a{x,p) ='^{xp + px) - {x){p) (2.8) 

evaluated in the initial state. In particular, for long times, the phase space spreading behaves 
according to 

^"^ (2.9) 

This means that any initial state becomes spread out in phase space on the (typically very 
short) timescale {H/'ykTy^'^. (See Ref. {25! for similar calculations). 

In the case above, a free particle, the spreading continues indefinitely. However, for a 
bound system with dissipation, equilibrium is eventually reached. For a simple harmonic 
oscillator at thermal equilibrium, for example, the ratio of thermal to quantum fluctuations 
is 

^ (2.10) 

h fvjj 

At typical laboratory temperatures and frequencies that are not unrealistically fast (for 
macroscopic systems), this number can be very large, of order 10^° say. These estimates will 
be relevant later on. In brief, they show that the density operator becomes very spread out 
in phase space (and hence slowly varying) very readily. 

This much is known material and will be used later. However, we now note that the exis- 
tence of the Wigner representation suggests an alternative method for constructing commut- 
ing position and momentum operators. It is well-known that when the Wigner function is 
sufficiently spread out, it becomes positive and may be regarded as a probability distribution 
for the variables p and q (which clearly commute, since they are numbers, not operators). 
The variables p and q do not correspond exactly to the operators p and x, but they are close 
when the Wigner function is spread out. In fact, it is easy to see that we have the following 
correspondences: 

qW{p,q) -{xp + px) (2-11) 

pW{p,q) ^ ^{Pp + pP) (2.12) 

That is, multiplication by g or p in the Wigner representation corresponds to the operation 
of anticommutation with x or p on the density operator (and it is easy to show that these 
two operations on p commute). The point here is that these operators on p are operations 
on the space of density operators which have no counterpart in terms of operations on pure 
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states. This is therefore not a route to producing the desired pair of commuting operators 
envisage by von Neumann. 



III. AN EXTENSION OF THE BALIAN-LOW THEOREM 

We first show that there is no projection operator satisfying the three properties 
Eqs. (11.181) . fll.lQp and fll.20p . This is an almost trivial extension of the proof of the Balian- 



Low theorem given by Battle 17 1. 



We consider the object Tt{Exp), which we assume exists. (If it does not, i.e., is infinite, 
this means that E has infinite dispersion in either x or p, so is not phase space localized). 
We have, from the three properties of Emn, Eqs. (ll.l7|l - (ll.l9p . 

Ty{Exp) = Tr (ExEmnp) 

mn 

= Y,^T [ExU^nEUl^p] 

mn 

= ^ Tr {EUran{x + na)E{p + hm)Ul^) 

mn 

= Tr {E_rri-n{x + na)E{p + hm)) 

mn 

= "^Tt {E_rn,~nXEp) 
mn 

= Tr {xEp) 

= Tt{Epx) (3.1) 

Or in other words, 

Tt{E[x,p]) = (3.2) 

which, via the commutation relations, implies that TrE = 0. Since E > 0, this means that 
E = 0. So no projector satisfying the three properties exists, if one insists that they be 
satisfied exactly. Hence, as indicated we will relax the conditions in what follows. 

We note in passing that this simple result has implications for the Balian-Low theorem in 
the pure state case. It may seem that one could avoid the Balian-Low theorem by dropping 
the requirement Eq. fll.lOp and requiring that the states iV'nm) are generated from more than 
one fiducial state. However, one could then use those fiducial states (assuming they are 
orthogonal) to construct a projection operator E satisfying the properties Eqs. (II. 170 - 01. 19p 
and the above result shows that the Balian-Low theorem is not in fact avoided. 
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IV. AN ALMOST COMPLETE SET OF ORTHONORMAL PHASE SPACE LO- 
CALIZED STATES 



We now show how to construct an almost complete orthonormal set of phase space local- 
ized states, from which we can construct the projector E. Our starting point is the set of 
states considered by Low [l^, 

^nm{x) = a-h{x - na) e^™/'^ (4.1) 

where h{x) a window function on [—la, la]. These clearly satisfy the three conditions, 
Eqs.f ll.lOj) . (11.131) and (I1.14p . The Fourier transform of these wave function is 

i^nmip) - (1) (p^/;^ _ 27rm) " ^^'^^ 

from which it is easy to see that (Ap)^ diverges because ipnmip) goes to zero like 1/p for 
large p which is not fast enough. Differently put, the window function in Eg. (14. II) causes 
the derivative of the wave function to involve the 5-functions 5{x — na ± |a), which are not 
square-integrable. These properties are fully in line with the Balian-Low theorem. 

However, inspired by a suggestion of Zak [26], it is easy to see from Eg. (14.21) that we 
may take linear combinations of these states which fall off like 1/p^ and therefore have finite 
(Ap)^. In particular, the states 

i^nmi^) = ^ {'^n,2m{x) + iJn,2m+l{x)) (4.3) 

have this property. One can also see this in configuration space: they vanish at x = na± ^a, 
and the offending (5-function appearing in their derivative therefore causes no problems. 
They are not complete, consisting of just "half" the states in the momentum direction, and 
indeed the states missed out are the states 

Xnlii.^) = {'^n,2mix) - 1pn,2m+l{x)) (4.4) 

and these have infinite dispersion in p. 

But it is not hard to see that the "remainder" states Xnm decay like 1/p in momen- 
tum space, so a further process of "halving" is possible to produce more states with finite 
dispersion. That is, we define 

= ^ (xl^L(^) - X^X+iix)) (4.5) 
13 



with new remainder states 

Xii = ^ (xi'L(^) + xi'L+i(^)) (4.6) 

So the set of states il^nm, V'nm, Xnm is complete and orthonormal, but the states ipnm have 
infinite dispersion. 

We can continue in this way to define the sequence of states 

i^Lrn^ ^ = (Xn,2m+1 " Xn,2m 

and 

(i^+l) ^ ^ f (K) (K) \ 

A.nm y^Xn,2m+l An, 2m j V .°/ 

for X = 1, 2, 3, • • • . If we truncate the sequence at some finite value ol K, K — N, say, then 
the set of states ipnm-, i^nm ■ ■ ■ i^nm together with the remainder states Xnm are orthonormal 
and complete. We therefore have the completeness relation 

TV 

E E ^i^H^)^im^(2/)* + E xilH^)xi^H2/)* = ^(^ - y) (4.9) 

K=l n,m n,m 

(K) 

In some sense, "most" of the states, namely the ipnm , have finite dispersion, and "some" 
of them, namely the Xmn have infinite dispersion. In this way, as N increases, the infinite 
dispersion anticipated from the Balian-Low theorem is pushed into a progressively smaller 
set of states. (An interesting question is whether the limit N ^ oo may be taken in any 
meaningful or useful way, but we will not pursue that here). 

Since each state is a linear combination of the ipnm, one may also derive the following 
general formula, 

2^-1 

« = 2-^/2 J] cf V'„,2-m+, (4.10) 

j=Q 



where 



and 



We also have 



cf = (-1)^ if < J < 2^-^ - 1 (4.11) 



cf = -(-1)^ if 2^-^ < J < 2^ - 1 (4.12) 



2^-1 

xi^) = 2-^/^E (-l)'>n,2-m+, (4.13) 

j=0 
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This may also be written, 

Xi^^ = + (4.14) 

The states are not all obtained from a single fiducial state, since we have 

li^i"^) = Ui^W) (4.15) 

where 

C/W = U,^,K^ (4.16) 

and 

1^^) = l^o^o) (4.17) 

There are therefore N fiducial states for the set ijj^ , K — 1, - ■ ■ N. 

At some length, one can compute the averages and dispersions of p and x in the fiducial 
states |'0'^). One obtains 

(.).^H?(.--i) (4.18) 

{x)k = (4.19) 
,2 /27r^\' (22-f^-l) , , 

(Ax)^ = g (4.21) 

(Note that (p)^' is not zero. Because there is more than one fiducial state, it does not 
appear to be possible to shift the states so that (p)k — in the fiducial states without 
spoiling orthogonality.) 

The construction described above is concisely summarized as follows: For each set of 2^ 
lattice points in the momentum direction, there are 2-^ — 1 states with finite dispersion and 
just 1 state (the remainder state), with infinite dispersion. 

Consider the following simple example to illustrate the construction. Consider the 8 
lattice points m = 0, 1 ... 7 in the momentum direction, so = 3. Then the 8 states which 
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depend only on these points are 



,(1) 1 



^^^l = + (4.22) 

V'n.l = 2 ^^"'4 ~ ^'^.S - ^n,6 + V'n,?) (4.23) 

^r^,0 = ^ (V'n.O - ^n,l + V'n,2 " V'n.S " ^"^,4 + V'n,5 " ^n,6 + V'n,?) (4-24) 
Xn'o = ^ (^n,0 - + ^n,2 " V'n.S + V^n,4 " V'n.S + V^n,6 " ^nj) (4.25) 

There are 7 finite dispersion states and one infinite dispersion state. It is easy to see that 
they are orthonormal and form a complete set on these 8 lattice points. Note also that 
widths of the states increases with K, but this is in some sense offset by the fact that 
the states become progressively sparser. 

So far we have been working with a complete set of states. We now come to the specific 
form of the proposal to relax the requirement of using a complete set of states (in the 
construction of commuting operators P and X): we quite simply drop the remainder states 
Xnm, which have infinite dispersion, and use only the incomplete set , K = 1,2 - ■■ N, 
which have finite dispersion. That is, in each set of 2^ lattice points in the momentum 
direction, we use 2^ — 1 out of the 2^ states. 

It seems likely that this approximation will be valid for sufficiently large A^, for suitable 
density matrices. We will show in the next section that reasonable results can be obtained 
using the ipi^ alone. Here, we briefiy look at the remainder states Xmn to see under what 
conditions they may be dropped. Note first that the states ipi^{x) have the property 
that they vanish at the end-points of the intervals, x = na ± a/2 (as they must, so that 
their derivative does not have a 5-function). Not surprisingly therefore, the remainder 
states Xnm{x) become narrower and progressively more concentrated about the end-points 
as — + oo, as one can see from Eq. (l4.14p (since they are in some sense making up for the 
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fact that the V'nm'' (x) vanish at the end-points, and the whole set of states is complete). This 
suggests that, under suitably coarse grained conditions, the behaviour of theses states at 
single points will become insignificant. We will see this more explicitly below. 

We may now use the almost complete set of states to construct the desired projection 
operators Enm localized on phase space cells. We choose the phase space cells to have 2^ 
lattice points in the p-direction and 1 lattice point in the x-direction. We have found 2^ — 1 
states with finite dispersion in each of those cells. We therefore take the projector E to be 

E (4-26) 

K=l m=0 

It is a sum over all states depending only on the lattice points to 2^ — 1 (in the momentum 
direction). Importantly, this is an exact projection operator which localizes onto a region of 
phase space - it satisfies 

E^ = E (4.27) 

exactly. The projection operator for any other cell is easily obtained by unitary displacement, 
using steps of size 2^ lattice points in the momentum direction (and single steps in the x- 
direction): 

Enm = UjZ^E{ujZ>y (4.28) 

These projectors clearly satisfy the exclusivity condition, Eq. fll.l7p . exactly, but do not 
satisfy the exhaustivity condition Eg. (11. 181) since we have 

+ E \x'n2){xi!2\ = ^ (4.29) 

nm n,m 

But, as we have argued, we expect the x terms to be negligible under suitable conditions so 
the projectors Enm should be almost exhaustive 

Y^Enm^l (4.30) 

nm 

The approximate nature of this property may not in fact matter for many practical appli- 
cations. The phase space projector onto a large cell V in phase space is defined by 

Et= Enm (4.31) 

The projector onto the region outside V is then defined to be 

Ey = 1- Ey (4.32) 
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so they are trivially exhaustive, Er + Er = 1. The key point here is that the remainder 
states dropped in the construction of E have infinite dispersion so they do not naturally 
belong in the construction of a phase space projector for a finite region of phase space. 
We will need some further properties of E. We have 

TrE = 2^-1 (4.33) 

which means it projects onto a phase space region of size (2^ — l){27rh). The object E/TiE 
may be thought of as a density operator and we can compute averages and variances to see 
the properties of E. We have 

(P). - ^ (4.34) 
2^fi /jA'-i _ I) (4,35, 



and 



2 



i^P)l - E ^ (4-36) 

K=l 



2^+%2;i2 



?>o? 

where we show only the leading terms of large A^. We also have 



(4.37) 



{x)e = (4.38) 
(Aa:)| = ^ (4.39) 

Since {p)e 0, it is in fact useful to perform a simple translation in momentum and define 
a related projector E' with all the same properties as E except that Tt{pE') = 0. We will 
use this in what follows. 

The construction of an exact projector with the above properties is the main achievement 
of this section and will be used to construct the commuting X and P operators below. 

Finally, we make two minor remarks. First, note that position and momentum enter in 
the construction in very different ways. However, the constant a is arbitrary, so may be tuned 
to make the width of the projector E arbitrarily small in either the x or p direction. Also, 
the states Eq. fl4.ip and (14. 2 p are a Fourier transform pair, so we could easily interchange 
them and start with a set of states that have perfect localization in p, instead of x. It would 
of course be of interest to find a construction in which x and p entered on an equal footing. 
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Second, we note for comparison that Omnes has made extensive use of phase quasi- 
projectors of the form 

Pr= dpdq \p,q){p,q\ (4.40) 

n 

where \p,q) are phase space locahzed states, such as the coherent states |9[. These have 
proved very useful for discussing emergent classicahty in quantum theory. However, in 
contrast to the projectors constructed here, these are not exact projectors, since they obey 
the approximate relation 

^ Pr (4.41) 

(and so Pr and 1 — are only approximately exclusive). It would be of interest to revisit 
some of Omnes results using the exact projectors constructed here. 



V. VALIDITY OF APPROXIMATE COMPLETENESS 

Now we come to a crucial check of our approach, which is to determine the conditions 
under which working with an approximately complete set of states gives reasonable results. 
The density matrix p of the system satisfies Trp = 1. If the states ioi K = 1, - ■ ■ N 

are approximately complete, then we should have 

N 

K=l nm 

which is the same as 

J2^HEnmP)-0 (5.2) 

nm 

We check this. It is most useful to exploit the Wigner representation, Eqs. (]2.ip . (]2.2p together 
with the property Eq. (12.31) . so we have 

E E (V'i^Vi^i^^) = E E / ^p^i ^^""^(p^ ^) ^'^(p + ^ + (5-3) 

K=l nm K=l nm 

where Wp is the Wigner function of p and W'^^'^ is the Wigner function of Now the 

crucial step is to approximate the discrete sum over n, m with an integral over continuous 
variables, p = q = na and we have 

ETTT , I^^^Tih , /" dpdq , _ . . ^. 

Wp{p + m,q + na) ^ / ^K+i^^f^ Wp{p + p,q + q) (5.4) 

nm 

^ (5.5) 



2^+%/i 
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This approximation is valid as long as the Wigner function of p is slowly varying over phase 
space volumes of size 2^(27r/l). Since K runs up to A^, we require slow variation on a scale of 
size 2^{27rh). This is typically the case for density matrices that are sufficiently decohered, 
as we saw in Section 2. We now have 

N N ^ 

- (5-6) 

K=l nm K=l 

= 1 - ^ (5.7) 

Hence we do indeed get a result close to 1 as long as is sufficiently large. 

One other related result is worth recording here since it will be used in the next section. 
Any density operator satisfies the relation 

/ ^f^'(^'^)/'f^(^'^) = l (5-8) 

where U{q, k) is the unitary shift operator in phase space. This is easily proved using the 
Wigner representation above, since we have 

{x\U\q,k)pU{q,k)\y) = J dp e^^^'-y^ _ fc, ^ - g) (5.9) 

and integrating the right-hand side over k and q yields 27ihS{x — y). 

One would expect that, for sufficiently slowly varying Wigner functions, a discrete version 
of the result Eg. (15.81) would hold. So suppose instead of the operator U{q,k), we use the 
operators Un^ defined in Eg. (14. 161) . We then have 

dp e^^(^-^) W,{p -—m, - na) (5.10) 

Using the same approximation and same steps as in Eg. (l5.4p to do the sum over n,m, it 
follows that 

E ^Ul^S)'pUl^ ^ i (5.11) 



2^ 

nm 

and therefore 



K=l nm 
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VI. CONSTRUCTION OF COMMUTING POSITION AND MOMENTUM OP- 
ERATORS 



We now use the projection operator E'^^ to construct commuting position and momentum 
operators, as outlined in the introduction. They are, 

X Kn. (6.1) 

nm 

P ^Y^Prn Kn. (6-2) 

nm 

where Xn — na and Pm — m x 2^ {27rh/a). Clearly 

[X,P]=Q (6.3) 

as required. 

We need to determine whether these operators are close to the original canonical pair, 
x,p. To do this we need some measure of distance, || P —p \\. We will define this by 



II P-p ||^= Tr (^(P-p)>j (6.4) 

(and similarly for X) where p is the density operator of the system, which we will assume 
is reasonably spread out in phase space. A more general measure of distance would include 
some sort of optimization over p, but, for the class of reasonably decohered density operators, 
the result depends only weakly on p. 
We have 

II P - p |P= Tr ((P2 - 2Pp + f)p^ + Tr (P\p, p]) (6.5) 

where the operators P are moved to the left in the first term in preparation for inserting the 
explicit form for P below (since 7^ 0). It is useful to introduce the notation tr(- • • ) to 

denote a trace over the incomplete set of states j'^im'')- We may then write 

II P - p f= tr (^(P2 _ + p2^p^ + 4 + Tr {P\p, p]) (6.6) 

where 

dl = Tr(pV) - tr(pV) (6.7) 

In the previous section we saw that Trp and trp are both very close (to order 2~^) and it is 
reasonable to expect the same if p is replaced with p^p, so dj^ will negligible and we drop it. 
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We also expect the term Tr (j^[p, p]j to be small, since a decohered p will be approxi- 
mately diagonal in momentum. Inserting the explicit expression for P, it is easily shown 
that 

Tr [P[P, p]) = J2 ^-Tr {[p, iu!,Z^)^ pU^Z^]) (6.8) 

nm 

Since we are working in the approximation in which p is sufficiently slowly varying that the 
sum over n becomes an integral. It is then easily seen (using Eq.( ]5.8t for example) that, in 
this approximation, the object 

T.(UiT)'pUnZ' (6.9) 

n 

is diagonal in momentum. Therefore the term Eq. fl6.8p is zero to the approximation we are 
using. 

Consider the one remaining term in Eg. (16.61) . Inserting the explicit for for P, we have 

\\P-P\\' = Y.^T{E'^^{Pm-pfp) (6.10) 

nm 

= Y.'r{muZ'ypuiz') (6.11) 



From Eg. (15. Ill) , we have 



and so 



E(^™^)W^^ - ^ (6-12) 



P-pf ^ (Ap)| (6.13) 



3a2 

where only the leading term for large is given. 
A similar (and simpler) calculation shows that 



(6.14) 



X - X 11^ ^ {AxYe (6-15) 

.2 



a 
12 



(6.16) 



Putting all these results together we obtain 

\\ P-p\\ .\\ X -x\\ ^ Ch (6.17) 

where 

r = 9^/2 

3\/2 
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C = 2^/2^ (6.18) 



This result is valid in the approximation of large N and for slowly varying density operators. 

Eqs. (15.71) and (16.171) show that needs to be chosen to be "large" to get approximate 
completeness, yet "small" for the commuting operators to be close to the original canonical 
pair. It seems likely, however, that there is a range of intermediate values that will meet both 
of these requirements. For example, take = 20. Then C ~ 10^, safely within the limit 
estimated in Eq. (ll.3p . so the difference between the old and new operators will be completely 
invisible to macroscropic observations. The error due to approximate completeness is of order 
, which is about 10~^. So there appears to be a large regime in which the approach works 
well. 

VII. PROBABILITIES FOR POSITION AND MOMENTUM 

The results Eq. (15.71) and Eq. (l6.17p are only indicative, and the true test of closeness of 
the old and new operators is a comparison of the probabilities for x,p and X, P. 

Consider the probability that the variable X lies in the range Ax, where Ax is the 
interval [nia,n2a], for some integers ni,n2. The probability is 

p(Ax) = Tr(pi?A,) (7.1) 

where the projector E/^^ is defined by 

^Ax= E JL^nm (7-2) 

(and the loose notation n G Ax means n G [^i,?t-2])- The probability for lying outside the 
region Ax is defined using the projector 1 — -Ea^j so that the probabilities add to 1 exactly, 
and the approximate completeness discussed earlier poses no problems. The probability for 
P is similarly defined, in terms of a projector E^^ defined by 

mgAp n 

Both of these probabilities therefore involve the probability pnm associated with our basic 
phase space cell of size 2^(27r/i), given by 

p.^ = Tr(K™P) = Tr {E' (UiZy)^ pU^^^) (7.4) 



23 



It is more usefully written in the Wigner representation, 

Pmn = 27r^ J dpdq We' ip, q) Wpip + P^r^, q + X„) (7.5) 

where, recall, X„ = na and = m x 2^ [27rh/a). 
The probability for the variable X is 

Since the Wigner function of p is assumed slowly varying, the sum over m may be approxi- 
mated by an integral over P„j regarded as a continuous variable. The p in p + P^ is absorbed 
into the integration and we obtain, 

p{Ax) ~ ^ ^ / dpdq WE'ip, q) p{q + ^n, ? + ^n) (7.7) 

Now the p integral may be carried out with the result 

p{Ax) ^aJ2 fdq P{q + ^n, g + ^n) (7.8) 

Since E' is phase space localized, one can see that the first part of the integrand is a smearing 
function peaked about g = and with width (Ag)^ a^/12. It is normalized to 1 (for large 
N) when integrated over q since TxE' = 2^ — 1. If we assume that 

A| » ^ (7.9) 

then the presence of the smearing function makes no difference and we obtain 

p(Ax) ~ a ^ p{na,na) (7-10) 

neAx 

If, as we assume, the density operator is sufficiently slowly varying for the discrete sum to 
become an integral, we obtain 

p{Ax) ^ [ dX p{X,X) (7.11) 



X 



the usual probability for position. 

Similarly for P, with analogous approximations, we obtain 
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^ E P^Pm.Pn.) (7.12) 

where p{p,p') is the Fourier transform of the density matrix p{x, y). Again, for slowly varying 
density operators the discrete sum becomes an integral and we obtain 

p(Ap) ^ / dPp{P,P) (7.13) 

JAp 

which coincides with the usual probability for p. 

We therefore find that the probabilities for X,P coincide with those for x,p as long as 
the following conditions hold: 

(i) is sufficiently large that the errors 1/2^ are tolerably small 

(ii) The density operator p is slowly varying on scales of size 2^(27r/l). 

(iii) The widths of the projections satisfy A^Ap ^ 2^(27r/i). 

The key restriction is (ii), the restriction on the density operator. Eqs. fl2.9p and fl2.10p 
indicate that the density operator can easily become sufficiently broad for (ii) to be satisfied. 
For example, for a time-evolving state, (ii) will be satisfied for 

*»(^) 2"/^ (7.14) 

This can easily be extremely short, even for large N. Similarly, for a state close to thermal 
equilibrium, (ii) is satisfied for 



^»2^ (7.15) 
nw 



which is again easily satisfied. 



VIII. SUMMARY AND DISCUSSION 

We have constructed a pair of commuting operators X, P which, at sufficiently coarse 
grained scales, are close (in a variety of ways) to the canonical position and momentum 
operators x,p. These commuting operators offer a new way of defining the relationship 
between approximate and exact decoherence. 

There are two ways in which this programme could be developed and improved. First of 
all, this paper has concentrated on the technicalities of constructing X and P. It would be 
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useful to develop more details of the conceptual framework in which they are used to discuss 
emergent classicahty. 

Secondly, the present approach works at the coarse grained scales of about 2^ phase 
space cells. Although arguably small to classical eyes, there are some ways in which this 
scale is quite large, and there are indications that some version of the present approach 
should work at finer scales. For example, it is known that the Wigner function (and also 
the P-function) become positive when coarse grained over just one or two phase space cells 
(rather than 2^ cells, as here), a feature that is often taken as indicator of approximate 
decoherence and emergent classicality |27,]- This suggests that there might be another way 
to construct commuting position and momentum operators which does not require such a 
large amount of phase space coarse graining. It would, for example, be of interest to see if 
von Neumann's original suggestion (involving an explicit orthogonalization of the coherent 
states) can actually be made to work. 

These and related questions will be pursued in future publications. 
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